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Abstract. Predictions of J/ifj formation at RHIC via "off-diagonal" combinations of 
charm and anticharm quarks in a region of color deconfincmcnt arc confronted with 
initial data from the PHENIX collaboration. We find that the measured centrality 
behavior places significant constraints on the various parameters which control model 
calculations of J /%p formation. Within present statistical and systematic uncertainties, 
one can map out a region of parameter space within which the contribution of formation 
in a deconfined phase is allowed. As these uncertainties decrease and new data from d- 
Au interactions becomes available, it is expected that definitive tests for the presence 
of this formation mechanism will be possible. We anticipate that the rapidity and 
transverse momentum spectra will prove decisive for a final determination. 



1. Introduction 

Initial data at RHIC energy on J/ip production in Au-Au collisions has been eagerly 
awaited, in terms of a signal for the presence of color deconfinement pQ. Of special 
interest is the possibility that a direct extrapolation of anomalous suppression from the 
SPS energy range could be supplanted by a new formation mechanism fueled by the 
presence of multiple pairs of charm quarks in each nuclear collision at sufficiently high 
energy [2] . One can argue on general grounds that the resulting J/ ip formation will grow 
quadratically with the total number of unbound charm quark pairs, with a normalization 
factor which depends on the specific formation mechanism jHj. Predictions for this type 
of formation have been made for two specific models. The statistical hadronization 
model j3] assumes that at the time of hadronization, charm quarks are distributed into 
hadrons according to statistical rules, incorporating an additional fugacity factor to 
conserve charm. The kinetic formation model E] considers J / ip formation within the 
region of deconfinement, and calculates the net number remaining at hadronization due 
to a competition between formation and breakup reactions. In this work, we confront the 
kinetic formation model predictions with the initial data from the PHENIX collaboration 
at RHIC. This initial data E] consists of a measurement of J/ip produced at central 
rapidity in Au-Au collisions at yfs = 200 GeV, presented in three centrality bins. In 
addition, there is data on J/ip production in pp collisions in three rapidity bins. There 
is also an indirect measurement of initial charm production via detection of high-p t 
electrons jH]. 
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In this presentation we first give in the next section a brief summary of the kinetic 
model formulation, with a discussion of primary uncertainties and ranges of parameters. 
The following section is devoted to a comparison of the predictions of this model with 
the initial PHENIX numbers, both the absolute values and the centrality dependence. 
Finally, we explore the regions of parameter space which are consistent with the central 
values of the initial data, and comment on aspects for the future. 

2. Kinetic formation model 

The kinetic formation model is most easily motivated in a scenario where suppression 
of J/ip due to deconfined color is related to their breakup via collisions with gluons [1QJ. 
Given the distribution of gluons and the cross section for this reaction, one can calculate 
the rate for destruction of J/ip in the region of deconfinement. It is then clear that one 
should also take into account the inverse reaction, in which charm and anticharm quarks 
interact to form J/ ip with emission of a gluon. The net effect is a competition between 
these two processes, which one can express via a Boltzmann equation for the charm 
quark and J ftp populations. 

= A F iV cPc - - A D iV J/v , Pg , (1) 

where p denotes number density, and the reactivity A is the reaction rate {av Te \) averaged 
over the momentum distribution of the initial participants, i.e. c and c for A^ and J/ijj 
and g for A^. In our calculations, this equation is solved numerically while enforcing 
exact charm conservation. At SPS energy, one finds that the formation process is 
negligble, since the average number of charm quark pairs produced even in central 
collisions is much smaller than unity. At high energy where one expects the number of 
initial charm pairs N c5 will become large, one can find an approximate analytic solution 
for Eq. [T] which exhibits the anticipated quadratic dependence. 

Nj/^t) = e(t) x [Nj/4t ) + Nl f A F [V(t') e^')]" 1 dt'), (2) 

J to 

where e(t) = e DPs would be the equivalent suppression factor in this scenario 
if the formation mechanism were neglected. Here we have assumed the densities are 
uniform within a deconfinement volume V(t). We use a thermal gluon density, initial 
temperature T as a variable parameter, and fixed final deconfinement temperature. 
The deconfinement volume is assumed to expand isentropically, for scenarios in which 
transverse expansion is controlled by a parameter vtr. The initial population of J/ijj 
is taken to be a fraction x of the initial number of charm quark pairs. The cross 
section is taken from an OPE-based model, based on the color dipole interaction of a 
nonrelativistic bound state with a coulomb bound state spectrum. For the centrality 
dependence, we utilize a calculation of the number of participants as a function of 
impact parameter. We also calculate the participant density in the transverse plane, 
and define an effective transverse the ratio of total to density. The participant 
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density is also used to specify the centrality dependence of initial temperature T . All 
centrality-dependent quantities are then scaled to their b=0 values. 
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Figure 1. Predictions for J/tjj formation. All use N c£ = 10 and x = 0.01 . 

Our initial calculations were performed for To = 0.3, 0.4, and 0.5 GeV, x = 0.01, 
and N c5 (b=0) = 10. To minimize the dependence on N cS , the ratio Nj/^ / N c£ 
was presented. Shown in Figure ^ are predictions for this ratio as published in Ref. 

where the legends are defined. The primary sensitivity was due to the charm 
quark momentum distribution. As might be expected, a thermal distribution was 
most efficient in formation, while distributions with various increasing rapidity widths 
predicted smaller numbers. All of the curves rise with increasing centrality, and have 
magnitudes near or above the assumed initial value (x = 0.01). 

3. PHENIX data vs. predictions 

For an initial comparison of these predictions with the PHENIX data, one must translate 
the measured dNj/^/dy at y=0 to the total Nj/^, and calculate the initial N c5 in each 
centrality bin. For the former, we assumed a flat J/i/j rapidity distribution with an 
effective Ay = 4. The N c5 for each centrality region were scaled to the N c5 (b=0) = 10 
assumed in the calculated curves, and varied with centrality according to the nuclear 
overlap function Taa(Jo) evaluated at an impact parameter corresponding to the central 
values of the experimental bins. The calculated ratios obtained in this way from the 
data are shown in Figure El where the error bars are taken from a sum of statistical and 
systematic uncertainties. One sees that the central values decrease with centrality. (This 
is the same behavior as shown in Reference |7] for the measured J/ip rapidity density 
scaled by the number of binary collisions.) It must be remembered that the ratios 
extracted from the experimental data in this way are only valid if the experimental 
value of total charm production proves to be the same as that assumed in these model 
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calculations. For ease of comparison, one set of model predictions for T = 0.4 GeV and 
charm quark momentum widths Ay = 1,2,3,4 are shown from among those in Figure 
[D Also shown are three additional model curves which use the exact charm quark 
momentum distribution from a LO pQCD calculation. 
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Figure 2. Indirect comparison of formation model with PHENIX data. 

These latter results are in general agreement with the preliminary data. For a 
more systematic comparison of the data with formation model results, we will restrict 
the charm momentum distributions to the pQCD set. 



4. PHENIX data and formation parameters 

In this section we adapt the model calculations to compare directly with dNj/^/dy 
data. The calculated Nj/^ have one component originating from initial production and 
a second from the combined formation/dissociation processes. The first is converted to a 
rapidity density according to the measured J / ip rapidity distribution from production in 
pp interactions [7] , and the second from a model calculation of the rapidity distribution 
which follows from the formation mechanism [11]. One also needs the N c5 values, which 
we estimate from preliminary open charm cross section 9\. Since there is presently an 
uncertainty of order factor of two in this estimate, we choose to present calculations for 
a wide range of iV c5 (b=0) = 5,10,15,20. We start with zero for our initial production 
parameter x. Given these parameter constraints, we scan over a range of T Q to select 
sets which are roughly compatible with the data. We show in Figure 01 such a set. 
One may question how the formation model is able to produce these curves, which 
rise with centrality less rapidly than binary scaling in order to be compatible with the 
data. The answer is that the formation model involves an inverse volume, which would 
normally provide a linear N part dependence. However this volume is time dependent 
and is integrated over the deconfinement lifetime along with factors sensitive to the 
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dissociation probability. We also extend the parameter space to include a nonzero 
transverse expansion of the deconfinement volume, shown in Figure 0] One sees that 
there is a significant correlation between sets of allowed parameters, and that a precise 
experimental constraint on N c5 , for example, would place significant constraints on the 
others. 




Figure 3. Kinetic model variation with Tq and N cS . 
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Figure 4. Kinetic model results for dNj/^/dj with transverse expansion. 

Finally, we show for completeness an alternate scenario in Figure EJ Along with the 
data and one set of formation model parameters which reproduces the central values 
of the measurements, we show two solutions in which only the dissociation mechanism 
is nonzero, using T = 0.3 GeV and 0.4 GeV which control the gluon density and the 
deconfinement lifetime. For comparison, the centrality dependence which would result 
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from pure binary collision scaling are also shown. It appears that even if the uncertain- 
ties in the data points were reduced considerably, one could probably find an acceptable 
dissociation fit. In that case, it will be essential to compare not only the magnitudes 




Figure 5. Comparison of formation formation calculation and binary scaling with 
PHENIX data. 



and centrality dependence, but also to confront the different scenarios with measured 
J/i/j rapidity and transverse momentum behavior. Work along these lines is underway. 
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